Delay Induced Excitability 
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We analyse the stochastic dynamics of a bistable system under the influence of time-delayed 
feedback. Assuming an asymmetric potential, we show the existence of a regime in which the 
systems dynamic displays excitability by calculating the relevant residence time distributions and 
correlation times. Experimentally we then observe this behaviour in the polarization dynamics of 
a vertical cavity surface emitting laser with opto-electronic feedback. Extending these observations 
to two-dimensional systems with dispersive coupling we finally show numerically that delay induced 
excitability can lead to the appearance of propagating wave-fronts and spirals. 
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The influence of time delayed feedback on the dynam- 
ics of otherwise Markovian processes with noise is cur- 
rently a very active topic of research, as memory effects 
have been found to be important in many systems in ar- 
eas ranging from biology to medicine and electronics . 
Even though many of the general properties of systems 
with memory are still unexplored, various methods for 
engineering of the delay to induce instabilities or stabi- 
lize complex dynamics have already been suggested Q • 

As a generic model to study the effects of delay on 
a bi-stable stochastic system one can consider a parti- 
cle trapped in a double well potential. In the absence 
of feedback the particle makes random transitions be- 
tween the two minima of the potential and two consec- 
utive jumps can be considered as being independed of 
each other. When feedback is present, however, correla- 
tions are introduced into the system and the dynamics 
becomes inherently non-Markovian. Let us consider the 
following standard stochastic differential equation 



dU 
dx 



ex(t 



(1) 



which describes an over-damped motion in a double well 
potential U(x) under the influence of memory. The 
length of the delay interval is r and e is the feedback 
strength. Here rj(t) describes Gaussian white noise with 
( v (t)} = and (rj(t)r](t')) = 2D6(t - f). In the following 
we will denote the delayed state as x T = x(t — r). 

Without delay (e = 0), the random switching between 
the two minima of the potential, xi, r , can be very well 
described by Kramers' theory [2j. The residence time 
distribution (RTD) Pi. r {T) in each well is then given by 



Pl,r{T) = 7i,r 



(2) 



where the inverse of the switching rate is the so-called 
Kramers' time, T^ r = 1/7;. r . It describes the average 
time between two transitions 



T, 



2n 



Here Xq is the position of the local maximum between 
the two minima x^ r and A£/j >r = U(xq) — U{xi >r ) is the 
potential barrier for each well. 

The statistical properties of eq. for e / have 
recently been analysed in the case of the generic sym- 
metric potential U(x) — —x 2 /2 + x 4 /4. Noting that 
the potential can be rewritten to include the delay term, 
V(x) — U (x) + ex x T , it can immediately be seen that the 
potential barrier, and therefore the escape rates, become 
dependent on the state at the earlier time t— r 0]. Hence, 
the continuous system of eq. can be approximated by 
a discrete two-state model, r(t) = ±1, with transition 
rates that depend on the delayed state []|. These rates 
depend on the relative value of r(t) and r(t — r) and are 
defined as p\ if r(t) — r(t — r) and P2 if r(t) 7^ r(t — r) 
for a symmetric potential. Tsimring and Pikovsky have 
recently derived the power spectrum for this model and 
shown the existence of coherence resonance |5jL Other 
works have calculated and measured the RTD 0, P 



Here we show that the presence of an asymmetry in 
the potential can lead to the appearance of excitability 
[S| • In an excitable system a perturbation above a certain 
threshold transfers an otherwise stable steady-state (the 
resting state) into an excited state (the firing state). The 
system then goes through a well defined refractory cycle, 
before returning to its initial state. During the refractory 




y/U"(xi, r )U"(x ) 



exp [AUi, r /D] . 



(3) 



FIG. 1: (a) Schematic definition of Kramers' residence times 
for an asymmetric potential without (upper graph) and with 
(lower graphs) negative feedback, (b) Schematic of the po- 
tentials during one excitable cycle. 
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cycle it is impervious to any external perturbation and 
cannot 'fire' a second time. Systems following this kind of 
dynamics are widespread in biology, chemistry and many 
other areas and there are several mathematical models 
existing to describe such dynamics • It is worth point- 
ing out, that for noise driven excitable systems different 
manifestations of a higher degree of order at a finite noise 
level have been discovered, for example coherence reso- 
nance ,9] and wavefronts and even self-ordered patterns 
like spirals in spatially extended systems • 

For the asymmetric potential the presence of feedback 
leads to a dependence of the amplitude of each potential 
barrier on the position of the particle at the time t— r and 
the dynamics has to be described by four characteristic 
times Tj^ r . The times TJ + and Tf~ are the mean escape 
times from the left well when the particle was in the left 
or right well at t — r, respectively, and T~ and T, + are 
the mean escape times from the right well if the particle 
was in the left or right well at t — t, respectively (see 

Fig.nt). 

Even though all transitions are inherently stochastic, 
we show below that for negative feedback (e < 0) a pa- 
rameter region exists, in which the systems displays ex- 
citable behaviour. To illustrate this let us consider the 
case of an asymmetric potential with T r <C Ti (for e = 0) 
such that the minimum corresponding to x < exists for 
any value of the feedback, x T , but the minimum for x > 
only exists for x T < (see Fig. Qd). Let us also assume 
that the particle has been in the left well for t < and 
it is driven by noise into the right well at t = 0. For 
small noise amplitudes {T~ 3> t) the particle will almost 
certainly remain in the right well until the value of the 
feedback term changes at t = r. The potential minimum 
on the right hand side then disappears and the particle 
will move toward the remaining minimum at x < 0. Af- 
ter a further time interval r, the potential will switch to 
its initial bi-stable shape and the above process can start 
over again. Such a cycle is characteristic of an excitable 
system. The initial state is stable under small pertur- 
bations but large perturbations can induce a jump. The 
firing corresponds to the transition of the particle into 
the right well and the refractory time consists of the two 
subsequent periods of duration r before the initial state 
is established again. During the first one the refractory 
time the particle resides in the well x > 0, from which 
transitions are unlikely and during the second one only 
one minimum is present and transitions cannot occur. 

To show that the dynamics described above can indeed 
be observed from eq. Q we have numerically calculated 
the RTDs in the excitable regime (see Fig. EJ. The po- 
tential was chosen to be 



U(x) = ~x 2 + ^x 4 + ^kx\ 



(4) 
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FIG. 2: Numerical RTDs for the left (full line) and right well 
(dashed line). The delay interval is r = 1. The inset shows 
the change in slope at T = 2 for particles in the left well, 
which is due to a higher order effect in the feedback strength. 



tern has fired) at the time t — 0, one can clearly see 
from Fig. [21 that the transition probability back into the 
left well (dashed line) is effectively zero in the interval 
[0, t] The spike appearing for T > r then indicates the 
non-stochastic roll-over to the left well once the feedback 
has changed and the potential possesses only the single 
minimum on the left hand side. From there no further 
transition is possible until at a time r later the potential 
changes back to its bistable shape (full line in Fig. |5J). 
After that the distribution for t > r follows Kramers' 
law. Fig. clearly shows the existence of a refractory 
time that forbids firing twice in an interval of 2r. 

Let us compare our continuous system to a discrete 
model. Although the simplest discrete model to describe 
excitability can be constructed using a two-state system, 
it was recently shown, that networks of excitable units 
display a wider spectrum of dynamical effects when they 
are approximated by three state models [lllfl2l |. In these 
models only the transition from the resting into the fir- 
ing state is stochastic, while the other two transitions 
are deterministic with fixed waiting times. It is clear 
from the discussion above that such a model is very close 
to our bistable system with feedback. Even though all 
transitions in the potential are at all times inherently 
stochastic, the two intervals during the cycle in which 
the transition probability is very low can be well approx- 
imated by a waiting time of length T w — t (see Fig. |5J|. 

To compare the power spectra we will assume that the 
stochastic output of the discrete model is described by a 
function that has a value 1 if the system is in the firing 
state and is zero otherwise. The spectrum can then be 
calculated analytically 0] 



2 — 2 cos(wr) 



(7 + 27 2 r)(I - cos(2wr) + £ sin(2wr) + 



2r 



with k = 0.3 and e = 1/3. Assuming that the parti- 
cle has made a transition into the right well (the sys- 



r) 
(5) 

Fig. [3] shows the power spectrum given by eq. J5J (up- 
per graph) and calculated from numerical simulations of 
eq. JU with the potential of eq. (@J (lower graph). Both 
curves are similar in shape and show their maxima close 
to odd multiples of ir and minima at even multiples of 
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FIG. 3: Power spectrum of the discrete (upper) and the con- 
tinuous model (lower). Here r = 1, e = 0.33, k = 0.3. 
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FIG. 4: Left: autocorrelation function of the system with 
increasing noise. For the uppermost curve D = 0.01 which 
increases downwards in steps of AD = 0.005. Right: correla- 
tion times. The maximum is a clear sign of ordering. 



7T. However, it must be noted that for the discrete model 
the coherence is an increasing function with noise [Tllll2| . 
while for the continuous model it is maximized for a fi- 
nite strength of the noise. To demonstrate this resonance 
we have calculated the autocorrelation function and the 
associated correlation times, tq, as a function of noise 
strength. The correlation time is defined as the integral 
over the square of the autocorrelation function || 



(x(t)x(t + f)) 



n 2 



(x(t) 2 



dt' 



(6) 



The left hand side of Fig. 0] shows the autocorrelation 
functions with noise increasing from the top to the bot- 
tom graph. It is clearly visible that the correlations are 
strongest for a finite noise level, which is confirmed in 
the plot on the right hand side, where the tc show a 
clear maximum for a noise amplitude of D ~ 0.04. This 
behaviour is closely related to the appearance of a phase 
transition from excitable to oscillatory media 



EXPERIMENTAL RESULTS 
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Excitability has been observed in the intensity of semi- 
conductor lasers with optical feedback 0. Here we ex- 
perimentally analyse the behaviour of the polarisation 
dynamics of a Vertical Cavity Surface Emitting Laser 
(VCSEL) with opto-electronic feedback as described in 
0, 0. The light emitted by a VCSEL is usually lin- 
early polarised but, as the injection current increases, 
the emission axis will switch between the two orthogo- 
nal states. Around this switching point one can find a 
current range in which the light polarisation randomly 
switches between the two axes. Previous experiments 



have been designed using these lasers to study Kramers' 
law 0], stochastic resonance and the properties of 
noisy time delay dynamical systems 0, 0] • It is worth- 
while noticing that the shape of the potential depends 
on the bias current applied to the laser. At the center 
of the bistable region both polarisations occur with the 
same probability and the associated potential is symmet- 
ric. The potential becomes asymmetric as the injection 
current is varied toward the boundary of the the bistable 
region and our laser was operated in this region in order 
to observe delay-induced excitability. 

Although spontaneous polarisation switching was ob- 
served without the addition of external noise, a noise 
generator was added in order to observe coherence res- 
onance. Fig. [S] (left) shows the measured RTDs for the 
lowest noise in the excitable regime where P x , P y refer to 
the orthogonal linear polarisation states. The important 
feature is the absence of switching events for times less 
than t. Additionally, both the stochastic nature of the 
firing transition and deterministic nature of the refrac- 
tory cycle can be seen in the spread of switching times 
being much less for P y than P x . Coherence resonance for 
this operating point is shown in Fig. (right) and quali- 
tatively reproduces the predicted numerical behaviour. 



TWO-DIMENSIONAL SYSTEMS 

The ability to sustain spatio-temporal patterns in sys- 
tems of coupled stochastic units is an important prop- 
erty of excitable systems 01- It was shown that wave- 
propagation and spiral patterns 0] are common in ex- 
citable systems and here we confirm their existence for 
spatially coupled elements as described by eq. We 
consider a situations where the elements are coupled via 
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FIG. 5: Left: Experimental RTDs for the x and y polarization 
of the VCSEL. Excitability is clearly visible from absence of 
residence times shorter than r. Right: Correlation times of 
the VCSEL output signal for increasing noise levels. 




FIG. 6: Left: Two dimensional wave-surfaces. The delay time 
in the left panel is twice the value of the delay time in the 
right panel. Right: Spiral pattern emerging from the initial 
condition (left panel) as described in the text. Black corre- 
sponds to the system in the firing state, grey to the refractory 
state and white to the normal state. The parameters here are 
r = 80, e = 0.343, k = 0.07 and D = 0. 



a diffusion term 

£)s 

— = s - ns 2 - s 3 - es T + V2L>£ + V 2 s . (7) 

Let us choose each element of the system to be initially 
in the resting state. At some point the noise will trigger 
one of the elements to fire and the coupling to neigh- 
boring elements will lead to short time correlations: in 
order to reduce the 'kinetic energy' between two neigh- 
bouring elements they will also undergo a transition into 
the firing state. After the refractory time the elements 
will transfer back in the initial state, marking the end 
of the wave-front. However, due to higher order effects 
in the feedback strength Q , the states at the end of the 
refractory time do not have the same Kramers' time as 
the initial ones. In fact, for negative feedback the new 
Kramers' time is much smaller and these elements are 
more prone to being excited again. This leads to the 
continuous creation of wave surfaces with wavelength r. 
In Fig. [Fj] (left) we show two examples of such wave sur- 
faces, where the delay time of the system on the right 
is twice as long as in the system shown on the left hand 
side. The two dominating colors represent the two re- 
fractory states for s > and s < 0. It is clear that these 
ordered waves can only be sustained at an optimal noise 
level. Too little noise will not lead to any excitations 
whereas too much noise will lead to fragmentation of the 
wavefronts. 

Instead of using noise to trigger the spatial order- 
ing, one can also choose appropriate initial conditions. 
Here we demonstrate the appearance of spiral patterns 
by preparing a system with one area where the elements 
have just fired after having been in the resting state for 
[— r, 0[ . A neighboring area has just made the transition 
back into the resting state, while having been in the re- 
fractory state for [— r, 0[. The rest of the system is in the 
resting state and has been there for [— r, 0] (see Fig. 
(right)). Integrating from these intial conditions leads to 
the evolution of a spiral as shown on the right hand side 
in Fig. E| 

In summary we have shown that including delay into 
bistable systems can lead to excitable behaviour. We 



have theoretically investigated the model of an asymmet- 
ric double well potentials and found clear signatures of 
excitability in the RTD. We have also compared the con- 
tinuous to a discrete model for excitability and found 
good agreement of the power spectra. We have calcu- 
lated the coherence times and shown the existence of 
coherence resonance. Experimentally we have demon- 
strated the above behaviour using a VCSEL in the bi- 
stable regime. For spatially extended coupled systems we 
have demonstrated the appearance of spatio-temporal or- 
der in the form of waves and spirals. 

This work was supported by Science Foundation Ire- 
land under contract number sfi/01/fi/co. 
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